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Abstract: What can be said about the spin of new particles without knowing their 
mass during an initial discovery phase at CERN? We consider this question in a topology 
where mass measurement is particularly difficult, pp — t- YY — t- IXIX, and introduce two 
new variables cos^^ and Afi which we prove are both independent of the mass of X. 
The variable cos^^ approximates the polar production angle of Y, and we find that it 
possesses greater statistical power in determining the spin of this particle than a previous, 
related variable, cos 9ii [1]. A^i is an asymmetry which can provide information about the 
couplings of a spin half Y. Because these variables can be used from the outset, without 
any knowledge of the masses of the new particles, we find here that it is possible to reverse 
the usual 'mass before spin' determination timeline. 
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1 Introduction 

To discover physics beyond the standard model at the LHC would be a great success. If 
new particles are created at CERN, their masses, spins and other quantum numbers will 
illuminate the details of the theory describing nature at the TeV scale. 

In this paper, we study a method of spin determination in the topology shown in 
figure 1: a new particle-antiparticle pair YY is produced, followed by the decays Y ^ IX 
and Y ^ IX where / is either an electron or a muon, and X and X are the lightest of 
the new particle states, which we assume are stable and invisible to the detector. We 
investigate the use of observables which do not take as input the mass spectrum of Y and 
X and show that the two introduced here are formally independent of the mass of X. The 
first is related to Barr's cos On [1], which approximates the polar production angle, cos^, of 
Y with respect to the beam axis in the YY centre of mass frame. We find an improvement 
on cos^ii's statistical power in determining the spin of Y. The second is an asymmetry 
which picks up on a term which can be present in the production matrix element squared 
for a spin half Y, |A^P cx ±cos6. We will discuss the possibility of using this asymmetry 
to provide information about the couplings of Y, something which is difficult to do using 
polar angle distributions alone. 

The problems that come with having two invisible particles in the final state of an event 
are well documented. Because mass measurement is usually 'easier' than spin measurement, 
in the sense that fewer events are required to obtain a reasonably accurate result, methods 
of spin determination are often happy to make the assumption that the mass spectrum 
of the theory is known. We do not go into detail here as to the methods which can be 
employed at hadron colliders to measure mass, and instead refer the reader to a recent 
review [2]. It is true, however, that the event topology we are considering is an example in 
which mass determination is particularly difficult [3]. This is down to not having two or 
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Figure 1. The event topology pp — > YY — > IXl X. 



more visible particles in the decay chain of y - there are no edges or endpoints which could 
come from the invariant mass distributions of such particles to provide relations between 
the new state masses and help to pin them down [2] . Nor is there an experimentally stark 
Mt2 'kink' [4-6] which would point directly to the mass of both particles. It would require 
many events to determine even the mass scale of the new physics, or else rely on mass 
measurements from elsewhere. At any rate, it is interesting to ask the question: what can 
we say about spin without any knowledge of the particle masses? 

Likelihood analysis techniques for spin determination use the complete event informa- 
tion and so can make statistically optimal statements, but their implementation becomes 
extremely complicated if the masses are unknown and need to be scanned over. Methods 
based on the total cross section [7] can also be very powerful, as can those which use kine- 
matic reconstruction [8-10], but they too require at least the masses to be known before 
they become useful. We shall investigate spin determination during an initial discovery 
phase, before these two techniques become viable, simply by studying the angular dis- 
tributions of visible particles in this topology - these are affected by the particular spin 
configuration at hand [1, 11-13]. We shall find that the properties of the observables intro- 
duced here make them ideal for use in a 'no mass input' approach and this leaves us in a 
position where it actually becomes possible to determine the spin of a new particle before 
we know its mass. 

In section 2 we discuss how the spin of Y can affect its production, identifying features 
which we can try to pick up on through the lepton momenta - this will also serve to 
introduce the set-up used in this paper. In section 3 we define the variable cos 6^, which is 
similar to cos 6u, and the asymmetry A^f^, which can be present if Y is spin half, and which 
gives insight into the vector-axial nature of its coupling. The interesting mass independence 
properties of both variables are studied here too. We demonstrate their use in determining 
the spin of Y in section 4 and present conclusions in section 5. 
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Figure 2. Production of YY by a spin one particle, Z, in the s-channel. 
2 Pair-production of YY 

For the present, we shall pretend that Y and Y are produced, do not decay and that we 
can measure their four-momentum perfectly. The following section will deal with the real- 
world problem of not being able to make this measurement when the final state cannot be 
reconstructed. 

We have chosen in this paper to illustrate the use of our variables by considering the 
s-channel production of Y and y by a spin one particle Z, as shown in figure 2, where Y is 
either a spin zero or a spin half particle. The particle Z could be either the standard model 
Z^ weak boson or it could be some new spin one state. We aim to discriminate between 
these two spin assignments where the Ys are produced according to distributions which 
go something like ~ sin^ 9 and ~ (1 + cos^ 9) respectively. We make the assumption that 
one can neglect all production diagrams other than the one in figure 2 - we shall discuss 
how this situation can arise shortly - and continue in a model independent fashion. We 
allow for the couplings between Z and Y to be those of the most general, renormalisable 
interaction: CZ^'{Y* {df,Y) - {df,Y*)Y) for a spin zero Y and Y^{V + A-f5)Y for a spin 
half Y, where C, A and V are the coupling constants. Because we are considering an early 
discovery phase at the LHC when the cross section for such a process won't be accurately 
known, we shall be comparing normalised distributions of variables. Given this, in the spin 
zero case there is no room for manoeuvre - the spin contribution is fixed as there is just an 
overall factor of the coupling constant, C. On the other hand, for the spin half case, there 
is the relative size of V and A, the vector-like and axial-like couplings of y to Z, which is 
a free parameter in this set-up. 

Under what conditions does the diagram of figure 2 dominate the production of YY7 
The importance of other diagrams, especially t-channel ones, must be minimal, since these 
can serve to make equal the angular distributions for spin zero and spin half Y [14] . Figure 2 
automatically depicts the dominant production mechanism in the case where Mz > 2My 
so that the particle Z, some new state, is resonantly produced. Alternatively, one could 
simply imagine that some symmetry of the theory describing these particles forbids any 
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t-channel diagram^. Because the former situation is more general, we choose to select some 
example masses for Y and Z so that Z can be resonantly produced, and use these values 
to illustrate the various steps along the way. We take My = 800 GeV and Mz = 2.5 TeV, 
and we shall return in section 4 to discuss how different choices for these masses affect 
our results (and will find that so long as Y and Y are reasonably well boosted in the 
LAB frame, which can happen for both Mz ^ 2My, then this is all we need to leave our 
conclusions unchanged). We shall also present results for the case when Z is the standard 
model Z". 

The two different spin assignments provide two different matrix elements to evaluate. 
We choose to do so in the centre of mass frame of the YY pair, where we write Y^s 
momentum as 



( 1 \ 



(i sin 6* 
\Pcos9 J 



(2.1) 



where /? is the velocity of Y and 7 = l/yl — /S^. In terms of these variables, the matrix 
element squared, summed over final particle spin and averaged over initial quark spin and 
colour is 



|7W|^ =AAF+ P(7) SPIN(/3,cose), (2.2) 

where 

^^^^ " (472 - M|/M2)2 + (TzMz/M^)^ ^^-^^ 

shows the effect of the Z propagator with a Breit-Weigner width Tz, N" is a numerical 
constant and -F-f is given in terms of the vector and axial couplings of Z to the quarks 

F+ = {f^f + Ut? ■ (2.4) 
Now we turn to the important difference. For a spin zero Y 

SPIN(/3, cos e) = (1 - cos^ 9) , (2.5) 

whereas for a spin half Y 

SPIN(/3, cos 6) = V^{2 + p\cos^ - I)) + (l + cos^ 6) ± 8-^:^VA /3 cos 9 .(2.6) 

We see an interesting final term in equation (2.6), asymmetric in cos^, which arises from 
an e^'^^'^ term coming out of the Dirac matrix algebra. This term is positive if one takes 
9 = in the direction of the quark, and is negative if it is taken in the direction of the 



^For example in [1] Z — Z'^ and so Mz < 2My and there is no resonance. Here the Fs are sleptons and 
conservation of lepton and baryon number forbids any tree level t-channel diagrams. 



-4- 



-1.0 -0.8 -0.6 -0.4 -0.2 0.0 0.2 0.4 0.6 0.8 1.0 -1.0 -0.5 0.0 0.5 1.0 

Cos 9 Cos e 

Figure 3. Two plots useful for spin determination if we were able to reconstruct the final state and 
thus know the momentum of Y. The left hand diagram shows the normalised distribution of cos 0, 
for the case when Y is spin zero and the case when Y is spin half with a purely vector coupling 
{A = 0), a purely axial coupling {V — 0), and mixed couplings {V = A ^ 1 is chosen). In the right 
hand diagram, the V = A case is also projected along the xi — X2 direction, revealing the 'hidden' 
asymmetry. Here, darker areas mean more events and the normalisation is arbitrary. We use the 
example masses for Y and Z described in the text. 




antiquark. Non-zero vector and axial couplings are needed at both the Zqq and the ZYY 
vertices for this term to be present in the first place. At this stage then we discuss the 
coupling of Z to the quarks, through and f^. The overall factor of in equation (2.2) 
will not affect the normalised distributions we are interested in and so the only place these 
couplings matter is in determining the relative size of the asymmetric term. We investigate 
by considering the case where /^/^/-F+ ~ C'(l) and find it useful to just take the values 
of these couplings as we would were Z the standard model Z^. This way, when we later set 
the mass and width of Z to those of Z^, this particle couples as the true standard model 
particle to the quarks. 

It remains to integrate these expressions over phase space and convolute this parton 
level process with the parton distribution functions [15] for the proton. We do this for pp 
collisions at 14TeV, using the values of our example masses outlined above. The resulting 
normalised differential cross sections are shown in the left hand diagram of figure 3 for 
four different cases: a spin zero Y; a spin half Y with a purely vector coupling (that is, 
setting yl = in equation (2.6)); a spin half Y with purely axial couplings (setting V = 
in equation (2.6)); and a spin half Y with mixed couplings, where we choose for illustration 
V = A = 1. Note the lack of asymmetry even when the VAcosO term is present. 

The prospect of distinguishing between a spin zero and a spin half Y looks promising 
using the variable cos 9 (remember we are imagining that we can measure Y and Y^s 
momenta and so know how to boost to their centre of mass frame), and we may even be 
hopeful about telling the difference between an axially coupling Y and the other spin half 
Ys. It is, however, disappointing that the asymmetry has been lost, especially as this could 
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give a handle on differentiating between the purely vector and the mixed axial- vector cases 
- something which is difficult to do using cos 6 alone. We should not be surprised that it 
is lost - the symmetry of the initial pp state of course prevents such an asymmetry being 
present in the distribution of cos^. To understand this in terms of equation (2.6), if we 
choose to measure cos 9 with respect to the direction of motion of the right-moving proton 
(which we call proton 1), then we are just as often defining = with respect to the 
direction of motion of a quark as we are an anti-quark. 

Although this asymmetry has been 'hidden', it is possible to measure it still - to do 
so we have to introduce an asymmetry of our own into the analysis. This can be done 
by distinguishing between events which are boosted in the positive z direction and those 
boosted in the negative z direction. To see how this can reintroduce the asymmetry, 
consider the case when the YY rest frame is boosted in the positive z direction in the LAB 
frame. This means that xi > X2, where xi and X2 are the momentum fractions of the 
partons coming from proton 1 and proton 2 (left-moving) respectively, and are given by 



XI = ^{Ey + Ey + p'y + p'y) 



X2 = ^{Ey + Ey - p'y - Py) , (2.7) 



1 

^' 

where ^/S is the centre of mass energy of the proton-proton system. For xi > X2, the parton 
distribution functions yield a higher probability for parton 1 to be a quark and parton 2 to 
be an antiquark, rather than the other way round. Overall there is a preference for the -|- 
sign in equation (2.6) and we see this in the right-hand diagram of figure 3 which plots the 
normalised cross section in the (cos 9, xi — X2) plane. We can now define an asymmetry as 

ATa::i>X2 , ajxi<X2 _ i\txi>X2 _ -\tX\<X2 
j^^^ _ ^'^cos6>>0 ^'^cos6><0 -'^cos6'<0 ^Vose>0 ^2 g) 

A^TOT 

where N^^^q^^ means the number of events with xi ^ X2 and cos 9 ^ and Atot is the 
total number of events. 

Both cos 9 and An can tell us something about the spin of Y and in the latter case, 
something about Ys coupling too. We shall now see whether we can construct similar 
variables using the lepton momenta which we could actually measure at the LHC. 



3 Letting the particles decay 

In this section we consider what happens when the pair-produced particles decay: Y — )■ IX 
and y — 7- IX. Because X and X are invisible to the detector, we are now in a position 
where we can no longer reconstruct the final state: cos0, xi, X2 and therefore An are no 
longer measurable. 

We let the above decay happen isotropically in the rest frame of Y , which means we 
take no spin correlations into account. To ignore spin correlations is to ignore an interesting 
way of determining spin in these situations (if the correlations are there, Y and Y are not 
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Figure 4. Comparing the distribution of the variable introduced here, cos 6*^ (left), with that of 
cosOii (right) for the different spin scenarios considered in the previous section. The faint, grey, 
dotted lines show the true cos 9 curves - those shown in figure 3 - for a spin zero Y, and for a 
spin half Y with a purely vector coupling (in the left hand plot, this lies beneath the purely axial 
coupling (green) curve). We use the example masses for Y and Z. The mass of X has absolutely 
no impact on the shape of these distributions. 



spin zero), but here they only have a small effect for two reasons. The first is that because 
all the final spin states are being summed over, the correlation is only observable through 
an overall production polarisation of Y and Y, if one exists. Secondly, Y and Y are 
produced with a sizeable boost in the LAB frame - this acts to wash out the effect of any 
lepton-lepton correlation. 

3.1 The visible cos^^ and the visible asymmetry Afi 

In reference [1], Barr introduces a variable which approximates cos 9 in the production and 
decay topology described above. This variable cos On = tanh(Ary;;/2) is invariant under 
longitudinal boosts as it is a function of Arju, the pseudorapidity difference between the 
leptons. This is an attractive quality at the LHC where we do not know the boost of the 
partonic centre of mass frame along the beam-axis. 

In a similar fashion, we now introduce a new visible variable cos^^, which is defined 
as the cosine of the angle between Vs direction of motion and the beam- axis in the centre 
of mass frame of the two visible leptons. That is one takes the sum of lepton momenta 
P = Pi + Pj, calculates the invariant mass rrin = ^/P^ and boosts the lepton momenta into 
the frame in which P = {mu, 0, 0, 0). This is the reference frame in which the angle 6^ is 
defined. This new variable cos 9^ coincides with cos 9ii when the leptons have equal and 
opposite transverse momenta. 

We plot the distribution of both of these visible variables in figure 4 for pp collisions 
at 14 TeV and using our example masses for Y and Z. Now we need to select a mass for X 
and we choose Mx = 500 GeV. The remarkable thing is that we could choose any value of 
Mx here {Mx < My, of course) and the distributions shown in figure 4 would stay exactly 
the same. We will discuss this rather interesting result below. In figure 4 the true cos 9 
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Figure 5. An observable version of the asymmetry shown in figure 3 (spin half Y, V = A = 1) 
is obtained by plotting events as a function of cos6'^ and — X2 ■ Darker areas mean more 
events and the normalisation is arbitrary. We use the example masses for Y and Z, and show the 
distribution for the choice Mx = 500 GeV. A larger value of Mx does squash this distribution 
towards xY — x^ = 0, but does not change the value of the asymmetry AYi- 



distributions are faintly shown and we can see that cos^^ tracks these true curves better 
than cos Oil, which suffers from a pseudorapidity 'drag down' at large |cos0//|. However, the 
important point is how well the visible curves can be distinguished from one another. By 
eye one might expect cos^^ to perform better - this is made into a statistical statement 
once detector effects have been taken into account in section 4. 

To observe the asymmetry we introduce two other variables xY and 

V 1 

X 



\ = —{Ei + Ej + pI+pI 



xY = ^{Ei + Ej-pt-pl) (3.1) 



and define the asymmetry 



^.V^^.V T^^T^ T^^T^ 

aV _ cose^>0"'" cos6"^<0 cose^<0 cos6''^>0 2) 

We plot this visible asymmetry in the {cos9^ ,xY — xY) plane in figure 5. 
3.2 Mass independence 

The variables we have just defined are exactly the same as the ones studied when we were 
pretending that we could measure Y and y's momenta, except that the roles of py and 
Py have been taken up by pi and p^. This is no surprise - it is down to the fact that both 
Y and Y have a relatively high boost in the LAB frame and so the leptons track their 
parents' direction well. We will begin by discussing this statement and will show that the 
angle between I and y's direction of motion in any frame of reference is independent of the 
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Figure 6. Plotting the relationship (3.5) between cosOy-i and cos9* (both defined in the text) 
for various values of /3'. A flat distribution in the centre of mass angle cos6'* leads to a distribution 
in cos 9y~i which becomes increasingly polarised towards cos 9y-i = 1 as 1" becomes more boosted. 



mass of X. The way in which cos^^ and are independent of Mx, a perhaps slightly 
more surprising result, follows easily from this observation. 

Consider the decay of Y in its own rest frame into X and /. 
momentum of / as 



Here we write the 



Pi 



E* 


E* sine* 
\E*cose* J 



(3.3) 



where E* = {My — Mj-)/2My. Now boost this to the frame in which Y travels with 
momentum pz = Myj'P' along the z axis (the primes are to distinguish these quantities 
from those used in the centre of mass frame of Y and Y) to obtain 



Pi 



(iE*{l + p'cose*)\ 


E* sin 61* 

V ys*(/3' + cosr) / 



(3.4) 



We can compute the cosine of the angle between / and Y's direction of motion 

7'(/3' + cosr) 



cos UY-l 



sin^ 0*+y2^/3' + cosr)2 



1/2 



(3.5) 



This result does not depend on our choice of the z-axis, and is true in a general frame of 
reference. There are two things we wish to point out. Firstly, as becomes large, the 
lepton really does start to track the direction of its parent's motion very well, as illustrated 
in figure 6. This is how the variables introduced above can pick up on the information 
about the polar production angle of Y . Secondly, E* has cancelled in this expression and 
so we see that this angle is independent of the mass of X. 
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Figure 7. Comparing the normalised distributions of cos 6^ (left) and cos du (right), after applying 
the kinematic cuts outlined in the text, for the four different spin scenarios considered throughout 
this paper. We use the example masses for Y and Z, and here X is not closely degenerate in mass 
with y - so long as this is the case, these curves are very insensitive to Mx (see text for details). 



The way in which E* cancels here points us to the reason why cos 6 and A[i are 
independent of Mx- The above result tells us that in any frame of reference, the momenta of 
I and I take the form pi = Ei{l,h) and pj = Ej{l, m), where n and m are unit vectors which 
have a distribution unaffected by the value of Mx- All dependence on Mx is contained 
within the overall factors Ei and Ej which are both proportional to AM^ = (My — M'^). It 
follows then that any dimensionless observable in any reference frame must be independent 
of the mass of X if the only dimensionfull inputs are components of the leptons' four 
momenta. The common factor of AM^ has to cancel. This completes the proof that cos 
is independent of Mx - We also point out that cos On , along with many other dimensionless 
variables which have been studied in the past such as A^^^, cos6'; — cos 6j, etc. all enjoy 
this property too. 

Turning to the asymmetry , we now have a dimensionfull input which is not a lepton 
momentum component: ^/S. We can see from equation (3.1) that xY — is proportional 
to the mass difference AM^, so that the distribution shown in figure 5 gets squashed 
towards = for larger Mx- However, the value of Mx has no effect on the sign 

of xY — xY , and so with the result for cos^^ we conclude that AY is also independent of 
Mx- 

These are nice features for variables if they are to be used without any mass input. In 
the next section we will see how they can be used to determine the spin of Y once some 
idealised detector effects have been taken into account. 

4 Collider simulation and differently boosted Ys 

We begin by declaring that this collider simulation will be one skewed towards the theo- 
retical end of the spectrum - in particular we consider only partonic processes and do not 
degrade the particles' momenta in any way. In reality many effects will lead to a smearing 
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Process 


CT(fb) 




WW 


0.98(1) 


-0.044(1) 


zz 


0.46(1) 


0.22(1) 


wz 


0.11(1) 


0.042(1) 



Table 1. Cross sections and values of AYi for the main standard model background processes. 
Exactly two opposite sign, same flavour leptons (e or /i) satisfying the cuts described in the text 
were required, along with missing energy. These were produced at leading order using the POWHEG 
Box [17] at the 14TeV LHC. Renormalisation and factorisation scales were set to the sum of the 
two vector boson masses in the process. The error shown is statistical only. In the ZZ case, only 
the two lepton, two neutrino final state was considered. 
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Figure 8. The cos6'^ distributions of the main standard model backgrounds in the two opposite 
sign, same flavour leptons plus missing energy channel, produced in the POWHEG box for the 
14TeV LHC with the cuts described in the text. 

of the distributions - for example, detector resolution and the issue of low pT initial state 
radiation^ - and we shall bear this in mind when making our conclusions. This approach 
is satisfactory because in [1] Barr carried out a realistic simulation far beyond the scope of 
this paper in the context of cos On and found that this variable could be used at the LHC 
for spin measurement ~ we can draw upon the similarities with this work below. 

We do however want to include detector effects which affect things at the parton level 
- this involves making kinematic cuts and studying the impact these have on the signal 
and the standard model backgrounds to this process. The cuts which we will employ 
are those of [1]: PT.i^ax > 40GeV, PT,i^in > 30GeV, where I 

max(min) is the Icpton with the 
larger (smaller) transverse momentum; a cut on the rapidity of the leptons, |7?z| < 2.5; and a 
requirement for missing transverse momentum pT,miss > 100 GeV. In addition, the following 
cuts are designed to remove the main standard model backgrounds: Mt2 > 100 GeV and 

^By this we mean immeasurable initial state radiation (ISR), in the sense that no hard jet is observed 
in the event. It would be interesting to further study whether hard ISR in the event could be turned into 
something useful, as was done in [16]. 
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Mil > 150 GeV'^. The main backgrounds for this channel of opposite sign, same flavour 
leptons are WW and ZZ production (both decaying to l'^l~'vu)^ and we will also consider 
the WZ production background. 

In figure 7 we plot normalised distributions for cos^^ and cos On, after these cuts 
have been made, for a spin zero Y and three types of spin half Y (vector coupling, axial 
coupling, and for couplings V = A = 1) at the 14TeV LHC and with our example masses. 
My = 800 GeV and Mz = 2.5 TeV. Because we are now cutting on the lepton momenta, 
a slight dependence on the value of the mass of X is introduced and so we quote that the 
distributions shown are with Mx = 500 GeV. However, for Mx in the range — 700 GeV 
these distributions only experience changes at the few percent level. For Mx above this 
range, the signal cross section becomes very small as the cuts really start to have an effect; 
in fact for Mx = 750 GeV the Mt2 cut entirely kills the signal^. If the case where X and 
Y are closely degenerate in mass is realised at the LHC, it will be very difficult to see the 
signal process above the standard model background anyway. The value for the asymmetry 

for a spin half Y with couplings V = A = 1 is found to be AYi = 0.16(1), where the 
error is purely statistical. 

The main standard model background cross sections can be found in table 1, along 
with values for the asymmetry, AY . Their distributions in cos 6^ - this time not normalised 
to unity - are plotted in figure 8. The importance of these backgrounds depends on the 
relative size of the beyond the standard model (BSM) process cross section. We include 
them to show their shape and magnitude. What was found in [1] was that they do not 
pose a major problem for using cos On to determine spin for BSM cross sections of order 
a few tens of femtobarns (provided X and Y are not closely degenerate in mass). We 
would expect this conclusion to hold true for cos^^ as well. Unfortunately, the standard 
model backgrounds, especially ZZ production, yield rather large values of A^. Since this 
has not been studied with realistic collider effects before, it is difficult to make a similar 
conclusion about the feasibility of observing and using this asymmetry at the LHC for spin 
determination, especially without a BSM cross section in mind. 

We now proceed to make a statistical statement about the power of cos in determin- 
ing the spin of Y, which comes down to being able to tell apart the different distributions 
in figure 7. Following the method of reference [19], we estimate the number of events N 
needed to disfavour one spin configuration S in favour of another, T, by a factor R as 

AT^J^g^ (41) 
KL(r,5)' ^ ' 

where KL(T, S) is the Kullback-Leibler distance [20] between the two curves. The results 
of this test are shown in table 2 for a value of R = 1000, where the number of events 
needed for distinguishing a spin zero Y from a spin half Y using cos 9^ are compared with 
those needed using cos On . We also show results for distinguishing an axially coupling spin 



^As the case we are considering involves an on-shell Z, an even higher cut on Mt2 could in fact be used 
to reduce the WW background further. 

*An animation which shows how these distributions change as Mx is varied from zero right up to the 
point at which the signal is killed by the cuts is available at [18]. 
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Variable Used 


T 


S 


cos 6 


cos 9^ 


cos Oil 


Vector 


Spin Zero 


28 


153 


179 


Axial 


Spin Zero 


25 


135 


158 


V = A 


Spin Zero 


21 


106 


124 


Vector 


Axial 


1126 


4008 


4870 


V = A 


Axial 


2487 


8863 


10794 



Table 2. This table shows the number of events N needed to discredit the spin configuration S in 
favour of spin T using the distributions of the variables cos 9^ and cos On in figure 7. For reference, 
the number of events needed using (the immeasurable) cos 9 are included. The top half of the 
table shows comparisons between spin zero and spin half cases, and the bottom half shows some 
comparisons between different coupling spin half cases. 
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Figure 9. This plot shows the number of events needed to discredit the spin zero Y hypotheses 
in favour of a spin half (purely vector coupling) F as a function of /^'^^ab ~ the average boost of Y 
in the LAB frame. This parameter is affected by the relative values of Mz and My, and is what 
is relevant for this method of spin determination (see text for details of how /S'^^ab varied here). 
Curves for cos0^, cos 9ii and (the immeasurable) cos 9 are compared. The results for the case when 
Z = Z'^, the standard model weak boson, are plotted as open squares. 



half Y from the other spin half cases (but not for telling apart a vector coupling and a 
mixed coupling Y - these distributions lie almost on top of each other in figure 7 and so 
will certainly be indistinguishable at the LHC). We shall comment on the number of events 
obtained in the context of spin before mass determination shortly, but for now we make two 
observations: firstly, cos^^ performs around 20% better than cos On in determining spin, 
and we would expect this to propagate through into a more realistic collider simulation; 
secondly, using cos 9^ to discriminate between the differently coupling spin half Ys requires 
enough events to render it an unrealistic approach at the LHC. Because of this, it would 
be very interesting if a more detailed and realistic detector study proved that Afi could be 
used to give information about the couplings of a spin half Y. 

We move on to address the question of how things change when we pick different values 
for the mass of Y and Z. The relevant parameter which is affected by these masses is the 
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boost of Y in the LAB frame, and in figure 9 we sliow how the number of events needed 
to disfavour a spin zero Y hypothesis in favour of a vector couphng spin half Y changes 
as the average value of this parameter, /3*^lab) varies. We choose to do this by keeping the 
mass of Y fixed at 800 GeV and varying the mass of Z. The lowest value of f3^l%^ is where 
we approach the limit of threshold production of YY in the Z rest frame. Shown also are 
results for the case when Z is the standard model Z^ - these points lie off the other curves 
because in this case the shape of the /3 spectrum is different, owing to the fact that we are 
no longer considering an on-shell resonance decaying to YY. As P^^X^ tends towards one, 
we see the behaviour we would expect from figure 6 ~ the number of events using the three 
different variables all tend towards the same value. This is because both cos and cos On 
approximate ever more closely the true cos 6 distribution, owing to the fact that the decay 
products from an increasingly relativistic Y will follow ever more closely the direction in 
which Y travels. On the other hand, as decreases, the number of events needed to 

determine y's spin using cos^^ and cos On rises exponentially, but we note that at the 
same time, cos 9^ becomes increasingly better at determining spin than cos On. In the case 
where Z = Z^, there is a factor of two gain in statistics in using cos 9^. The number of 
events needed diverge at threshold production due to the fact we considered an isotropic 
decay of Y. An animation which shows the changes taking place in the cos 9^ and cos 9ii 
distributions as /S'^lab is varied is also available at the website [18]. 

Is it possible to determine spin before mass here? Given what was said in the intro- 
duction about the measurement of mass, this question really reduces to whether it would 
be possible to determine spin in this topology during the course of an LHC run. Realistic 
collider effects will certainly increase the number of events indicated in table 2 and figure 9, 
but so long as Y and Y are produced with a boost greater than ~ 0.75 then one might 
expect to be able to make a spin measurement (spin zero or spin half) with a few hundred 
events. This would require the BSM cross section to be at least around the 10 fb level. 
Furthermore, X and Y can't be too closely degenerate in mass if the signal is not to be lost 
in the standard model background. These conclusions are backed up by what was found in 
the particular cases studied in [1]. Finally, with reference to early spin measurement, we 
point out that the value of R = 1000 chosen in equation (4.1) could be relaxed somewhat 
which would give a (less certain) indication of the spin of Y with fewer events. 

5 Discussion and outlook 

In this paper we have introduced two new variables which can be used for spin determina- 
tion at the LHC: cos^^ and . We considered the s-channel pair-production of particles 
Y and y by a spin one particle Z, and looked at the topology where YY — )• IXIX. We 
showed how cos 6^ can give information about the distribution of the polar production 
angle of Y, illustrating its use in discriminating between the case of a spin zero and a spin 
half Y. We found that in this theoretical study cos possesses more statistical power than 
the related variable cos 9ii, and would expect this to propagate through when more realistic 
collider effects are taken into account, as they were in the original study of cos 9ii [1]. The 
asymmetry which we introduced has the potential to identify a spin half Y which has 
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both vector and axial couplings to the Z production particle, something which cannot be 
done using cos0^ alone. Until further study, we reserve judgement on the usefulness of 
this variable for spin determination at the LHC - without specifying a particular value for 
the cross section of the new physics it is unclear how much of a role the standard model 
background would play in obscuring such a measurement, and a simulation of realistic col- 
lider effects is required here as well. It would be very interesting to see if a measurement 
of this asymmetry were possible at the LHC outside of the context of spin determination, 
especially as the standard model predicts a non-zero value for it in the two leptons plus 
missing energy channel. 

We proved that neither of these variables depend on the mass of the invisible decay 
product X, and pointed out that a number of other variables also enjoy this property - 
we are not aware that this was known before. The event topology is one in which mass 
determination is particularly difficult and so we considered the question of what could 
be said about spin without the knowledge of any new particle masses. Because of their 
independence of the mass of X, the question of whether these variables can be used to 
measure spin could be assessed on the basis of one parameter ~ the boost at which Y 
and Y are produced in the LAB frame. We found that so long as they are reasonably 
well boosted, discrimination between the two spin assignments in this set-up using cos 
becomes possible after a few hundred events, before any of the particle masses are known. 
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